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This  period  is  the  first  annum  under  this  new  grant.  Our  efforts  are  focused  on  the 
following  activities: 


1. 


Continue  to  study  the  8-operator  based  stability  test  algorithms,  due  to  the  importance 

f  this  subject.  We  have  performed  a  sensitivity  analysis,  and  have  shown  the 
ivantages  of  the  8-operator  based  stability  tests  developed  by  the  PI  under  the 
revious  grant  (N00014-90-J-1017).  Specifically,  we  have  defined  a  new  kind  of 
ensitivity  matrices  and  have  shown  that  the  8-operator  based  tests  have  sensitivity 
natrices  that  approach  those  of  the  continuous-time  tests  as  the  sampling  interval 
vanishes,  whereas  those  of  the  traditional  shift-operator  based  tests  grow  without 
bound  under  the  same  condition.  This  establishes  the  numerical  superiority  of  the  8- 
operator  based  stability  tests  analytically.  Submitted  papers  A1  and  A2  in  the  attached 

publications  report  are  summarizing  this  result.  Another  aspect  is  the  singular  case 

of  the  8 -operator  based  stability  tests.  Similar  to  the  shift-operator  based  stability 
tests,  the  8-operator  tests  developed  so  far  are  only  for  the  regular  case.  They  will  fail 
in  the  singular  case.  Therefore,  we  have  developed  a  stability  test  for  the  singular 

type  I  8-Schur-Cohn  test.  We  have  also  shown  that  the  limit  of  this  test  as  the 

sampling  interval  vanishes  is  a  corresponding  singular  test  in  the  continuous-time. 
Further,  we  have  reconciled  two  seemingly  different  approaches  in  dealing  with 
singularity  type  II,  one  in  the  shift-operator  discrete-time  and  the  other  in  the 
continuous-time,  through  the  8-operator  test  for  singularity  type  II.  Submitted  paper 
A3  is  for  this  result. 


2.  We  have  also  studied  the  normalized  lattice  structure  and  its  continuous-time  limit.  It 

turns  out  that,  unlike  any  other  filter  structure,  the  normalized  lattice  has  a 
continuous-time  limit  as  the  sampling  interval  vanishes.  This  explains  the  well- 
known  good  numerical  properties  of  the  normalized  lattice.  However,  this  limit  does 
not  realize  an  arbitrary  stable  transfer  function  in  the  continuous-time.  We  thus 

modified  the  normalized  lattice  so  its  limit  can  now  realize  an  arbitrary  stable  transfer 
function  in  the  continuous-time.  Various  aspects  of  stability  of  the  newly  modified 
normalized  lattice  are  studied.  The  results  are  summarized  in  A4  of  the  attached 
publication/presentations  report. 

3.  We  have  also  begun  to  investigate  the  t.^KD-based  recursive  least  squares  (QRD-kLs; 
algorithm  and  its  8-operator  version  for  adaptive  signal  processing.  We  have  obtained 
a  8-operator  based  algorithm  for  off-line  processing  (non-recursive)  and  its 
continuous-time  limit.  Currently  we  are  working  on  the  recursive  version.  We  are 
also  looking  into  finite  precision  implementation  aspect  of  the  QRD  and  8-QRD,  in 
order  to  understand  and  study  the  numerical  advantages  offered  by  these  algorithms. 

4.  Further  work  has  also  been  done  in  the  “generalized  delta  operator”,  in  cooperation 
with  the  Swedish  researchers  at  the  Systems  and  Control  Group,  Uppsala  University, 
Sweden.  Particularly,  we  are  interested  in  the  limiting  case  where  continuous-time 
autoregressive  model  parameters  are  to  be  estimated  from  the  discrete-time  data.  Our 
previous  work  in  using  a  modified  least  squares  (LS)  method  resulted  in  Publication  II 
of  the  attached  publication/presentations  report.  In  our  current  work  we  use  the 
conventional  LS  method  rather  than  the  modified  LS  method,  and  then  correct  the 
resulting  parameter  bias  at  the  last  stage  by  a  simple  scaling.  The  advantage  of  the 
new  approach  is  that,  since  the  conventional  LS  method  is  used,  we  can  now  use 
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recursive  algorithms  to  obtain  computational  savings.  We  have  submitted  a 
conference  paper  on  this  result  so  far.  But  more  publicatioil  is  on  the  way. 

5.  Blind  equalization.  This  work  was  supported  jointly  by  the  AASERT  program  under 
Grant  N00014-93-1-1032.  We  have  constructed  a  new  family  of  cost  functions  for  blind 
equalization.  Although  not  all  convex,  they  all  seem  to  be  unimodal.  Publication  12  in 
the  attached  publication/presentations  report  is  one  result.  More  work  has  been 
performed  since  then.  We  have  tried  to  show  the  unimodality  analytically,  which 
proved  to  be  very  difficult.  We  are  now  working  on  the  noise  issue,  and  are  also 
preparing  a  journal  paper  submission.  On  the  other  hand,  we  are  also  investigating 
the  potential  of  using  the  subspace  system  identification  methods  in  this  problem. 

The  results  are  very  preliminary  and  will  be  reported  subsequently. 

6.  During  this  period  we  also  revised  some  of  the  previously  submitted  papers.  This 
includes  a  finite  precision  paper  using  a  stochastic  differential  equation  approach 
(now  accepted  for  publication),  two  5-operator  stability  test  papers  (now  also  accepted 
for  publication),  and  a  paper  analyzing  conditioning  of  information  matrices  for  8- 
operator  based  algorithms.  All  these  papers  will  be  duly  reported  when  they  appear 
in  the  respective  journals. 
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